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ABSTRACT 



Aims. We aim to determine the constraints that observed chemical abundances put on the potential role of mass loss in metal poor 
dwarfs. 

Methods. Self-consistent stellar evolutionary models that include all the effects of atomic diffusion and radiative accelerations for 28 
chemical species were computed for stellar masses between 0.6 and 0.8 Mq. Models with an initial metallicity of Zq = 0.00017 and 
mass loss rates from lO^'^Moyr"' to 10^'^Moyr"' were calculated. They were then compared to previous models with mass loss, as 
well as to models with turbulent mixing. 

Results. For models with an initial metallicity of [Fe/H]o = -2.31, mass loss rates of about 10^'^Moyr"' lead to surface abundance 
profiles that are very similar to those obtained in models with turbulence. Both models have about the same level of agreement with 
observations of galactic-halo lithium abundances, as well as lithium and other elemental abundances from metal poor globular clusters 
such as NGC 6397. In this cluster, models with mass loss agree slightly better with subgiant observations of Li abundance than those 
with turbulence. Lower red giant branch stars instead favor the models with turbulence. Larger differences between models with mass 
loss and those with turbulence appear in the interior concentrations of metals. 

Conclusions. The relatively high mass loss rates required to reproduce plateau-like lithium abundances appear unlikely when com- 
pared to the solar mass-loss rate. However the presence of a chromosphere on these stars justifies further investigation of the mass-loss 
rates. 

Key words. Diffusion — stars: chemically peculiar — stars: mass-loss — stars: Pop II — stars: evolution — stars: abundances 



1. ASTROPHYSICAL CONTEXT 

■ Population II stars can be used as tracers for exploring the uni- 

■ verse's chemical evolution. Still, to correctly relate the chemical 
properties of the early Universe to the composition of the oldest 
observed stars, one must be able to map the evolution of their 
current abundances back to their original values. Doing so could 
help elucidate numerous questions pertaining to Big Bang nu- 

' cleosynthesis, the nature of the first supemovae, and the age of 

' the Universe. 

The nearly constant Li abundance of most low-metallicity 

] ([Z/H] <-1.5) halo field stars whose Tgs is betwe en 6300 K and 
5500 K - otherwise kn own as the Spite plateau dSpite & Spitd 
119821: ISpite et alj[l984l) - has puzzled astronomers for decades. 
Originally, many inferred that this plateau was representative of 
the Universe's original Li abundance, though it was promptl y 
shown that this could not be the case (Mich aud et alj [l984l) . 
Precise observations of the cosmic microwave backgrounqj 
have now led to a primordial lithium abund ancqj of A(^Li)- 
2.65+gg^, or even higher, A(^L i ) ^ 2 72+ggj (ICvburt et al., 2003; 
ISteigmanL 120071: ICvburt et all. feoOSi) . while the lithium content 
in metal poor dwarfs, ranging from A(^Li) = 2.0 to 2.4 is 



a fact or of 2 - 3 lower (iBonifacio et al.L 2007t Asplund et al. 



' From the Wilkinso n Microwave Anisotropic Background, WMAP, 
see lSpergel et al.ll2007i 
2 A(Li) = log[iV(Li)/A'(H) + 12]. 



2006: Charbonnel & Primas, 2005: M elendez & Rarmrezl 12004 
Bonifacia,2002:,Rvan et al...,1999.) . 

Changing the standard Big Bang nucleosynthesis 
model is a possible explanation, though it is hypothetical 
(Ichi kawa et al. 2004: Coc et al. 2 004: Cu mberbatch et al. 20071: 
Reg is & Clarksonll2012h . Even though Popula tion III stars may 
have destroyed some of the primordial lithium (" Piau et al.l2006l) . 
this cannot explain the large discrepancy. This leads many to 
assume that the solution stems from within the stars themselves. 
Possibl e stellar explanations may include destruction via gravity 
waves (iTalon & Charbonnell 12004 : Charbonnel & Talon 2005b, 
and rotational mixing ('Vauclair' 1988: Charbonnel et al. 1993 
Pinsonneault et al. 1999), as well as settling via atomic diffusion 
(Deliyannis et al. 1990: JProfitt & MichaudI 119911: ISalaris et al] 
.2000: .Richard et al...2005h . A higher Li content of subgiant (SG) 
stars compared to main-s equence (MS) stars in globular clus ters 
(GC) such as NGC 6397 (iKom et al.ll20"07l:lLind et al.l |2009l) fa- 
vors a model in which lithium is not destroyed, but simply sunk 
below the photosphere prior to dred ge-up. This result is, how- 
ever, sen sitive to the Tes scal e used (iGonzalez Hernandez et al.l 
2009: No rdlanderetal.1120121) . 

Additional constraints may stem from claimed observations 
of ^Li (Asplund et al. 2006), which, with a thermonuclear de- 
struction cross section 60 times larger than the heavier isotope, 
would severely constrain the stellar destruction of ^Li. These ob- 
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serva tions remain very uncertain (ICavrel et al.ltZOOTl iLind et all 
l2012h . Many recent studies have also claimed that below 
[Fe/H]^ -2.5, Li abundances lie significantly b elow the Spite 



self-co nsistently with radiative accelerations fr om [Richer et al 



1998) and corrections for redistribution from iGonzalez et al 



plateau abundance and with increased scatter (lAsplund et al 



20061 lAoki et all l2009l ISbordoneet al] l2010i ISbordone et a" 



20121) . Other observations suggest that there might instead be a 



lower plateau abundance for extremely metal poor (BMP) stars 
( Melend ez et al. 2010). The situation is complex and not cur- 
rently understood in detail. 

Although the surface convection zone (SCZ) of galactic 
Population II stars homogenizes abundances down to depths 
at which atomic diffusion is relatively slow, the long lifetimes 
of these stars, which reach the age of the Universe, allows for 
the effects of chemical separation to materialize at the surface. 
However, models that allow atomic diffusion to take place unim- 
peded lead to lithium und erabundances that a re too large at 
the hot end of the plateau jMichaud et alj[l984 ). By enforcing 
additional turbulent mixing below the SCZ, models of metal 
poor stars with gravitational settling and radiative accelerations 
are able to reproduce many observations of halo and globu- 
lar cluster star s (Richard et al., 2002b; VandenBerg et al., 2002; 
iRichard et all l2002at iRichard et al., 2005; Kom et al., 2006). 
These same models were als o successful in e xplaining abun- 
dance anomalies in AmFm dRicher et al.l 2000; Richard et al. 
1200 1!) and horizontal branch stars ("Mi chaud et al.ll2007l l2008). 

In Vauclair & Charbonnel (1995), gravitational settling is 
inhibited by unseparated mass loss at the constant rate of 
~ 3 X lO^'^Moyr '; this leads to near constant lithium abun- 
dances as a function of Teff, even at the hot end of the plateau. 
However, these calculations, which are discussed in Sect.|4] did 
not include the effects of radiative accelerations, nor did they 
follow atomic diffusion in detail for elements heavier than Li; 
theref ore, further inspection is warranted. 

In lVicketal1 (2010. 2011), hereafter Papers I and II, stel- 
lar evolution models with mass loss were introduced and shown 
to reproduce observed surface abundance anomalies for many 
AmFm and chemica lly anomalous HAeBe binary stars. In 
iMichaud et alj ( 1201 lb it was then shown that a mass loss rate 
compatible with the observed Sirius A mass loss rate leads to 
abundance anomalies similar to those observed on Sirius A. To 
further validate the model and to help constrain the relative im- 
portance of mass loss with respect to turbulent mixing in stars, 
the present mass loss models will be applied to Pop II stars and 
co mpared to models with turbulence from Richard et al. (2002b) 
and lRichard et al.l (l2005h . which are able to reproduce observa- 
tions of galactic halo stars, as well as of globular cluster stars on 
the subgiant and giant branches. 

In the following analysis and in our calculations, mass loss is 
considered in non rotating stars, since in these stars it might be 
the only process competing with atomic diffusion. To our knowl- 
edge, mass loss has never been observed in Population II MS 
stars; therefore, mass loss rates will be constrained solely via 
surface abundance anomalies. In Sect.|2] the evolutionary calcu- 
lations are presented, while the results are discussed in Sect. [3] 
In Sects. m and |5] the models with mass loss are compared to 
previous evolutionary models of metal poor dwarfs with mass 
loss, as well as to models with turbulence. In Sect.|6]results are 
compared to subgiant and giants both for models with turbulence 
and with mass loss. Conclusions follow in Sect.|7] 

2. Calculations 

The stellar evolution models were calculated as explained in 
Paper I and references therein. Atomic diffusion is computed 



[T995I) and iLeBlanc et al. (2000). Chemical transport is com- 
puted for 28 che mical speci es (all those included in the OPAL 
database. Jglesias & Rogersll996. plus Li, Be, and B). The mod- 
els are evolved from the pre -main-sequence (PMS ) with the 
abundance mix prescribed in iRichard et al.l (l2002bl) for Zq = 
0.00017 (or equally, [Fe/H]= -2.31). Most models were calcu- 
lated with an initial helium fraction o f 0.235 (corresponding t o 
the Yo used for calculations for M 92 in lVandenBerg et al.l2002l) . 
Since one of our major aims is to compare with similar mod- 
els calculated with turbulence, this value was chosen to facili- 
tate comparison (since it is the same as in IRichard et al.ll2002bl 
and lRichard et al.l2005h . Some models were also calculated with 
Yq = 0.2484, whic h was derived from WMAP obser vations 
( Cvburt et al.L 120031) and is di scussed in § 3.3.1 o f Richar d et al.l 
(|20()5'). As shown in Fig. 4 of lVandenBerg et al.l (119961) and Fig. 
5 of IRichard et al.. (2005.) . respectively, this small change in ini- 
tial helium fraction has little effect on the luminosity-versus-age 

relatiorH^ or on lithium isochrones. With respect to the solar 
abundance mix used in Papers I and II for Population I stars, 
the relative abundances of alpha elements are increased as ap- 
propriate for Population II stars (IVandenBerg et al.]l2000i) . The 
introduction of mass loss and its impact on transport is outlined 
in Paper I. 

The convection zone boundary is determined by the 
Schwarzchild c riterion with e nergy flux in the convection zone 
computed with iBohm- Vitensd (Il992h mixing length formalism. 
Abundances are assumed to be completely mixed within con- 
vection zones. No extra mixing is enforced outside of convec- 
tion zones. The effects of atomic diffusion materialize as soon 
as radiative zones appear as the models evolve toward the MS. 
These effects become more important as the radiative zone ex- 
pands toward the surface, where atomic diffusion becomes more 
efficient, and the convection zone becomes less massive, thereby 
reducing the dilution of anomalies. 

The unseparateqj mass loss rates considered range from 
10"'^ to lO^'^Moyr"'. Owing to uncertainties related to the 
nature of winds for Population II stars (see also discussion in 
Sect. 4.2 of Paper I), we have chosen to limit our investigation to 
unseparated winds in order to avoid introducing additional ad- 
justable parameters. 

In order to fully constrain the effects of m ass loss, we have 
chosen to use the same boundary condition dKrishna-Swamvl 
1966') and mixing length parameter, a - 2.096, as in Papers I and 
II (see also Se ct. 2 of Paper I andLV andenBerg et al. 2008). Most 
models in R ichard et all (l2002b(ai) and IRichard et al.l (120051) 
were computed with the Edd ington gray atmosph ere surface 
boundary condition (see also iTurcotte et al.l Il998l) . The main 
consequence is to change the value of a. In both cases, the mix- 
ing length is calibrated using the Sun so that very little difference 
remains in solar type models. Therefore, though the different 
boundary conditions might lead t o a different evolution of SCZ 
depths and masses (see Fig. 9 of IRichard et al.ll2002bh . the re- 
maining slight difference does not play a role in models with tur- 
bulence since the separation that modulates surface abundances 
usually occurs deeper As shown in Sect.|3] the same is true for 
the models with mass loss presented in this paper since separa- 
tion also occurs deeper than the bottom of the convection zone. 



^ In this context, unseparated refers to the abundances in the wind 
being the same as in the photosphere. 
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Fig. 1. The Heitzprung-Russel diagram, as well as the evolution of Teff, logg, the temperature at the base of the SCZ (Tsscz), the 
outer mass at the base of the SCZ (AMbscz), and the mass fraction of hydrogen in the core {X^) for stars of 0.6 - 0.8 M© with 
Zo = 0.00017 and M = lO^'^Moyr-^. 



3. Evolutionary models 

Evolutionary models of masses ranging from 0.6 to 0.8 M©, an 
interval of masses that spans the Spite plateau, are shown in a 
Hertzprung-Russel diagram (H-R) in Fig.[T] All the models have 
an initial metallicity of Zq - 0.00017, and a mass loss rate of 
10"'^Moyr"'. Figure[T]also shows the evolution of Teff, gravity, 
the depth (in mass, AMbscz, and in temperature, Tbscz) of the 
bottom of the SCZ, and the central H mass fraction. The 0.79 
and 0.77 Mq models are respectively at the MS turnoff at 12 Gyr 
and 13.5 Gyr, the assumed lower and upper limits of the age of 
halo dwarfs. 

For all models, the variations in surface convection zone 
depth have an impact on surface-abundance evolution since the 
efficiency of atomic diffusion varies with depth. Throughout the 



MS, AMbscz decreases until H is exhausted in the stellar core 
(i.e. until tumoff). Although this and other parameters depend 
on the mass loss rate, the effect is slight (see Fig.|2|. Just before 
the MS turnoff, there is a ~ 125 K difference in Teff between the 
0.8 Mo model with diffusion only and the 0.80W1E-12 modefl 
while the difference falls to about 15 K between the diffusion 
only and 0.80W1E-13 models. The temperature and mass at 
the bottom of the SCZ are related to both the He and Fe con- 
centrations in the SCZ. For instance, the temperature is lowest 
when both He and Fe have a small abundance at turn-off in the 
0.80W1E-13 model. The temperature is highest in the models 
with largest He and Fe abundances (and highest mass loss rate) 



0.80W1E-12 refers to a 0.80 M,, model with a mass loss rate of 



1 X lO^'^Moyr" 
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Fig. 2. The impact of mass loss on the evolution of the depth (in temperature and mass) of the SCZ in 0.8 Mq models with atomic 
diffusion only, and with mass loss rates from 10 '^Moyr ' to 10 '^Moyr '. upper row: HR diagram and the He and Fe surface 
abundances, lower row: Zoom on the Tetj, and on the mass and temperature at the bottom of the SCZ. 



and intermediate in the diffusion only model where the low He 
abundance is partly canceled by the high Fe abundance. The iron 
abundance at turn-off is largest (+0.47 dex) in the model with- 
out mass loss and decreases progressively to -0.2 dex as the mass 
loss rate increases to lO^'^'Moyr"' and then is very nearly nor- 
mal in the larger mass loss cases. The effect on stellar age at the 
turnoff, or equivalently, on the stellar lifetime, is less than the 
line thickness in the plots. While AMbscz/A/* varies by a factor 
of 100 (log AM/M» goes from ~ -2.3 to -4.3) on the MS for a 
0.80 Mq model (Fig.[TJ, the maximum effect of adding mass loss 
is below 0.08 dex. 



3.1. Radiative accelerations, mass loss rate, and internal 
abundance variations 

On very short time scales, convection homogenizes abundances 
from the bottom of the SCZ to the stellar surface; therefore, 
the competition between radiative acceleration (gmd) and gravity 
is only important below the surface convection zone. In Fig.|3] 
the internal variations of giad, velocity, and mass fraction are 
shown for Li, Mg, Ca, Ti, Fe, and Ni in the 0.80W1E-12 modefl 
Having gi-^^, velocity, and mass fraction aligned for a given ele- 
ment facilatates the discussion of the link between those quanti- 
ties. Figures for gy^^ and mass fraction of all elements are shown 
in App.lAl they allow to show the systematic variations of those 
quantities. The changes in g^^d over time are due principally to 
structure changes, since the abundances in Population II stars 



with [Fe/H]o = -2.31 are never large enough to cause satura- 
tioifl 

The internal variations of an element's gi-^d can be related 
to its ionization state as may be seen in the online App.lAl H- 
like electronic configurations lead to local maxima, while He- 
like configurations lead to local minima. As we move toward 
elements of higher atomic number, the H-like state is encoun- 
tered at greater depths. Carbon for instance, is in an H-like con- 
figuration at logAM/M» ~ -3, while for Si, this occurs near 
log AM/Mt ~ - 1 . Some elements, including Fe and Ni, never 
quite reach the H-like statqj. Unlike the H-like maximum, the 
second local giad maximum, which is clearly distinguishable for 
all elements heavier than S, is usually greater than or equal to 
gravity. This g^i increase corresponds to configurations between 
F-like and Li-like. 

For many elements, including He, Li, and O, the radiative ac- 
celeration is always smaller than gravity. For a mass loss rate of 
10"'^ Moyr"' , these elements are underabundant (bottom row of 
Fig-EJ throughout the star except close to the center Other ele- 
ments such as Be, Mg, and Fe have local g^fi maxima just below 
the SCZ, which, at 11 664 Myr, causes their surface abundances 
to be larger than just below the convection zone. All abundance 
variations caused by atomic diffusion are smaller than 0.2 dex 
(in contrast to the ones caused by thermonuclear reactions which 



^ This ma ss was chosen to facili tate comparison with the models and 
figures from lRichard et al.l(l2002bh . 



* The mass loss rate also has little to no effect on ^faj for this same 
reason, and because its effects on the structure are very small (see 
Fig.[3. 

' To be more precise, in the stellar core, these elements arrive at 
an equilibrium between multiple ionization states, including the hydro- 
genic one. 
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Fig. 3. upper row: Radiative accelerations (solid lines) and gravity (dotted lines); center row.- Wind velocities (dashed lines) and 
diffusion velocities toward the exterior (solid lines) and toward the interior (dotted lines); lower row: Interior concentrations as a 
fraction of original ones. These quantities are shown at four different ages for a 0.8 M© model with M - 10"'^ Moyr"' . The vertical 
lines indicate the position of the bottom of the SCZ. The ages are identified in the figure. 
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Fig. 4. The impact of mass loss rate on internal abundance profiles for 0.80 Mq models. A model without mass loss and models 
with five mass loss rates from 10"''*Moyr ' to 10 '"Moyr"' are shown and are identified in the figure. The internal concentration 
profiles are shown just before tumofi^near 11. 5 Gyr Atomic diffusion affects mainly the outer 10"^ of the stellar mass. Conventions 
are similar to those of Fig.|3] 



are encountered near and around the stellar core for Z < 8). For rapidly is the location of the bottom of the SCZ. The time varia- 
this mass loss rate, all elements are underabundant below the tion of interior concentrations is continuous and relatively small. 
SCZ, and down to about log AM/M» ~ -0.5. What varies most Larger variations are linked to local extrema of grad- For instance 
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X(Fe) has a local maximum at log AM/M, ~ -3 where grad(Fe) 
has a minimum and the diffusion velocity is inwards. These ef- 
fects will be discussed below in relation to the effect of varying 
the mass loss rate. Below log AM/M» ~ -0.5, elements heavier 
than O become slightly overabundant due to gravitational set- 
tling. Therefore, for these elements, there are underabundances 
of up to a factor of 1 .6 over the outer half of the stellar mass and 
overabundances of a few percent over the inner half. 

For models with lower mass loss rates, internal abundance 
variations become much greater (Fig.0. In the 0.80 M© model 
with a lO^'-'Moyr"' mass loss rate, Ca and Ti are locally over- 
abundant by a factor of up to ~ 3, while "^He, ^Li, Mg, and Fe 
become -0.8 dex underabundant. These variations can be linked 
directly to variations in grad- For example, the Ni underabun- 

dance that begins at logAM/M, 2.5, and reaches up to the 

BSCZ, is caused by the increase in its gi-^^ over the same depth 
interval. In this region, Ni is pushed into the SCZ, thus caus- 
ing the Ni abundance to be relatively larger at the surface and 
to decrease below the SCZ. In contrast, Ca is overabundant near 
log AM/Mt = -3.5 and above, since gmdiCa.) has a local maxi- 
mum just below this depth. 

In fact, when mass loss is added, one may distinguish two 
types of internal solutions: those dominated by atomic diffu- 
sion and those dominated by mass loss (see also Sect. 5.1.1 of 
Paper I). For any given model with mass loss, an element's in- 
ternal abundance solution is determined by the mass loss rate 
wherever the associated velocity, Vwind, is greater in amplitude 
than the settling velocity v'sett- For example, in Fig.|4] one first 
notes that the depth at which |vwindl - ksettl, which we call the 
point of separation, depends on the element and on the stellar 
age (through structure changes). For the 0.80W1E-12 model, 
mass loss dominates the internal solution of 'Li from the sur- 
face down to logAM/M» ~ -1.2, while this occurs nearer 
logAM/M, ~ -0.9 forNi. As discussed in Sect. 5.1.1 of Paper 
I, the consequence is that local abundances above the point of 
separation are not determined uniquely by the local g^^d, but in- 
stead adjust in order to conserve the flux arriving from just below 
the point of separation. If the wind-dominated solution extends 
down to AM = M x f, the surface abundance reflects giad/g at 
AM. If we come back to the previous example, at lOGyr — or 
lO'^yr — and for a mass loss rate of lO^'^Moyr"^ the surface 
reflects g^i^d/g at ~ 10"^^ M,. 

In general, a higher mass loss rate has the effect of re- 
ducing internal anomalies. For example, a mass loss rate of 
10^'^Moyr"' reduces the Li and Mg anomalies below the con- 
vection zone by ~0.6dex (Fig.|4]i. However, the abundance 
anomalies of, for instance, Fe do not always have a monotonous 
variation with mass loss rate. As one reduces the mass loss from 
10^'^Moyr"', the Fe abundance in and below the SCZ first de- 
creases and is 0.6 dex smaller below the SCZ in the 10 '"* Moyr ' 
model. In the SCZ it is 0.2 dex smaller However, in the model 
with a mass loss rate of lO^'^Moyr"', the surface Fe abundance 
is 0.35 dex larger than the original abundance, and without mass 
loss, it is 0.4 dex larger than the original one0. This surprising be- 
havior is related to the minimun of §iaci(Fe) at log AM/M, ~ -3. 
For a mass loss rate of lO^'^Moyr"', the wind carries away the 

surface Fe pushed by ^rad(Fe) above log AM/M* 3 since in 

one Gyr, it carries log AM/M» ~ -4 and the SCZ is slightly 

less than log AM/M» 4. The Fe flux that crosses the region 

where the Fe diffusion velocity is negative, at logAM/M, ~ 
-3, is reduced, leading to the observed underabundances. The 



lO^'^^Moyr"' wind carries a lower Fe mass from the outer re- 
gion and reduces the surface Fe by only 0.05 dex compared to 
the no-mass-loss model. 

On the other hand, Ca in the 0.80W1E-13 model in Fig.|4]has 
a larger surface overabundance than in the no-mass-loss model. 
This is because the Ca abundance must adjust to carry the flux 

strongly pushed upwards at log AM/M, 3 where gradiCa) is 

much larger than gravity. In the 10"''*Moyr"' model, the wind 
velocity is too weak to carry the flux when the diffusion veloc- 
ity is toward the interior, and there is a sharp peak of X(Ca) at 
logAM/M, = -3.5. The profile of the concentration is clearly 
linked to the profile of the downward oriented diffusion velocity, 
from Mg to Fe. 

3.2. Surface abundance variations 

For a given metallicity, surface abundances depend on stellar 
mass and age as well as on mass loss rate. The age and stel- 
lar mass dependance are demonstrated in Fig.|5] which shows 
the evolution of surface abundances for selected models with 
a mass loss rate of lO^'^Moyr"'. First, for this mass loss rate, 
all elements are underabundant throughout MS evolution, and 
until dredge up advects layers which are affected by thermonu- 
clear reactions. The underabundance amplitudes increase as stel- 
lar mass decreases, since settling becomes more important for 
models with longer lifetimes. Over a large fraction of the stellar 
lifetime, the overall shape of the surface abundance solution is 
quite similar for all masses since the overall shape of g^^d does 
not vary significantly over this mass interval. See Fig.|2]for the 
effect of varying the mass loss rate on the He and Fe abundances. 
At a given age, the mass loss rate affects the surface abun- 
dance of the various elements differently, since supported ele- 
ments are expelled through the surface, while sinking elements 
encounter advection by Vwind- For 0.80 M© models near tumoff 
(Fig.|6ll, a mass loss rate of lO^'^Moyr"' reduces He, Li, and O 
underabundances obtained in the diffusion only model by a lit- 
tle more than 0. 1 dex, while also reducing overabundances of CI 
and P by about 0.8 dex. It is then clear that for this mass loss rate, 
the wind's effect is greater for overabundant (i.e. supported) ele- 
ments, since the settling velocity still dominates for unsupported 
elements. For yet other elements the abundance variations are 
more complex: the overabundances of Fe and Ni become under- 
abundances, while Ca and Ti become more overabundant for a 
short period of time. For Ca and Ti, this is because Vwind has an 
effect just below the SCZ, where gi-^d - g (Fig.©; therefore, 
instead of accumulating below the SCZ, they are pushed into it. 
Similarly, mass loss advects the layers just below the SCZ, where 
Fe and Ni abundances are smallest, into the SCZ (see Fig.|4]and 
IA.3b . For the two higher mass loss rates, the surface anomalies 
are reduced to factors mostly below about 0.2 dex. While the 
0. 80W 1 E- 1 2 has generalized underabundances, the 0. 80W5E- 1 3 
model has small overabundances of elements from K to Mn (see 
Fig.lA3Tl. 



3.2.1. Surface lithium abundances 

Figure]?] shows lithium and iron isochrones for the diffusion only 
model, for models with mass loss rates of lO^'^Moyr"' and 
5 X lO^'^Moyr^', as well as for T6.09 modelfl at 12 Gyr and 



** We have verified that the 10 '^ M^yi ' model is the same as the 
no-mass-loss model. See Fig.]2] 



' The parameters specifying turbulent transport coefficients are in- 
dicated in the name assigned to the model. For instance, in the 
T6.09D400-3 model, the turbulent diffusion coefficient, Dj, is 400 
times greater than the He atomic diffusion coefficient at log Tq = 6.09 
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and varying as p ^, or 



Dj = 400£»He(ro) 



piTo) 



13.5 Gyr. To simplify the discussion only the effect of mass loss 
is considered i n this s ection. The T6.09 isochrones that are from 
Richard et al. ('2002b') are discussed in Sect.|5]and where the ef- 
fects of mass loss and turbulence are compared and the meaning 
of T6.09 is briefly described. 

In the diffusion only model, large Fe and Li anomalies are 
obtained around turnoff: lithium decreases by a factor of up to 
10 at 12 Gyr, while Fe becomes overabundant by a factor of 5. 
For the mass loss rate of 5 x lO^'^'Moyr"', the maximum lithium 
underabundance is reduced to -0.35 dex, whereas Fe overabun- 
dances are transformed into underabundances by factors below 
-0.25 dex. For the 10"'^ Moyr"' model, the maximum ampli- 
tudes are further reduced. While the anomalies in the diffusion 
only model near turnoff are much smaller at 13.5 Gyr than at 
12 Gyr, this effect is not as important for lithium, and is unno- 
ticeable for both models with mass loss. 

Below ~6000K, the lithium isochrones with and without 
mass loss differ only slightly for pretumoff stars — with lithium 
depletions between -0.12 and -0.14 dex near Tf^g = 5400 K — 
since the SCZ is too deep for mass loss to have an effect. 
Though all three isochrones have decreasing lithium abundances 
for models hotter than ~6000 K, the underabundance amplitude 
decreases as the mass loss rate increases. At 12 Gyr, the lithium 
abundance begins increasing again near 6200 K and 6400 K for 
the models with lO^'^Moyr"' and 5x10"'-' Moyr"' respectively. 
This is due to SCZ retraction, which exposes layers for which 
Vwind is dominant. 

The diffusion only lithium isochrones do not exhibit plateau- 
like behavior On the other hand, the mass loss isochrones are 
much flatter: the 5 x 10"'"^ Moyr"' model has a scatter of about 
~0.1dex and the 10~'^Moyr"' model has a scatter of about 
~0.05dex with respective plateau average values of about - 
0.3 dex and -0.2 dex below the initial ^Li abundance. 



(1) 



To simplify writing, T6.09 is also used instead of T6.09D400-3 since 
all models discussed in this paper have the D400-3 parametrization. 
The 0.80T6.09 model is a 0.80 Mq model with T6.09 turbulence. 
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4. Comparison to previous calculations with mass 
loss 

In previous evolutionary calculation s for Po pulation II stars with 
mass loss, 'V auclair & Charbonnell (1 19951) (VC95) found that 
both mass loss rates of 10"'^-^ (^ 3 x 10"'^) and lO^'^Moyr"' 
lead to lithium abundances that increase monotonously (by at 
most 0.5 dex for the lower mass loss rate) with T^ff for six mod- 
els ranging from 0.60 to 0.80 M©. They suggested that the higher 
mass loss rate completely wipes out the effects of atomic diffu- 
sion, and leads to a primordial lithium abundance of A(Li) = 
2.5 ±0.1, which they determined from the surface abundance of 
their hottest model. 

Unfortunately, a detailed comparison to their results is im- 
possible, since, even without mass loss, the stellar models are 
too different. These differences include Yq, Zq, ckmlt, the sur- 



face boundary condition, and the diffusion of metals. As a result, 
the SCZ evolution is very different. For instance, by comparing 
the internal lithium profile of the diffusion only model in Fig.|8] 
with the one from Fig. 4 of I Vauclair & Charbon nel (1995), one 
sees that at 10 Gyr, the SCZ (i.e. the flat part of the profile for 
both solutions) is about four times more massive in their mod- 
els. Furthermore, their 0.75 Mq models are still on the MS (with 
their Tetf still increasing) at 18 Gyr, while our 0.75 Mq models 
have already evolved toward the subgiant branch at 15 Gyr The 
difference in Yq is not large enough to explain this significant 
discrepancjo 



'" Though not shown, this was verified by evolving models with Yq = 
0.2484, whose lower H mass fraction in the core led to their passing 
tum off only 1 Gyr earlier than mod els with Yq = 0.2352. The models 
from lVauclair & Charbonnell (1 19951) have Ko = 0.2403 
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Fig. 8. [Left panel] Internal lithium profiles for 0.75 M© models with and without mass loss. Three models of different masses 
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IVauclair & Charbonnell(ll995l) . 



That being said, some general results can still be compared. 
For instance, at lOGyr, VC95 found that the surface lithium 
abundance for the lO^'^'^'Moyr"' model was 0.12 dex larger than 
in the lO^'^Moyr"' model, since the stronger mass loss rate 
reaches into the lithium burning layers. Our models also sug- 
gest that the strongest wind drags some material from the lithium 
depleted layers outward. Indeed, from Fig.|8] the wind veloc- 
ity for the 10^'^Moyr"' model at lOGyr dominates the settling 
velocity down to log AM/M, = -1, the point of separation, 
which is well into the burning layers. However, at lOGyr, this 
has not yet reached the surface since at lO'^yr. ^ mass loss rate 
of lO^'^Moyr"' has only exposed to the surface layers which 
are above ~ lO^^M,. In the same way, mass loss advected all 
mass above the point of separation, thus explaining the differ- 
ence in the depths at which Li falls to zero in the diffusion only 
and lO^'^Moyr"' models. 

Furthermore, the shape of the internal solutions in both sets 
of models is very diff'erent; for 10 '^Moyr^', VC95 obtain a 
nearly flat solution from the BSCZ down to Li-buming layers 
(see their Fig. 4), whereas our solution suggests an important 
atomic diffusion-induced gradient below the SCZ (see Fig.|8]l. 
In order to verify that the difference was not simply due to dif- 
ferences between SCZ masses during evolution, this result was 
compared to smaller mass models with the same mass loss rate 
(the 0.73W1E-12, 0.70W1E-12 and 0.68W1E-12 models). Even 
in models with more massive SCZs, the abundance gradient re- 
mains, since it is formed below log AM/M, ~ -2, the point 
above which the wind has imposed the conservation of flux. The 
reason for these differences between the two calculations is not 
evident. 

The Li profiles in the left panel of Fig.|8]are understood using 
the right panel and the assumption of unseparated mass loss from 



the surface. In the 0.75 M© diffusion only model, the gradient 
between log AM/M, = -1.8 and -2.7 is determined mainly by 
the increasing settling time scale as log AM/M, increases, that 
is closer to stellar center In the 0.75W1E-13 model the wind 
only shifts 10"^ of a solar mass in lO"* years which only very 
slightly modifies the profile around log AM/M, ~ -2. In the 
0.75W1E-12.5 model the efi^ect is slightly more pronounced as 
the mass loss velocity becomes a little larger, in magnitude, than 
the diffusion velocity as seen in the right panel. In the 0.75W1E- 
12 model however the mass loss velocity is larger than the dif- 
fusion velocity by a factor of up to six, and flux conservation 
determines the Li abundance between log AM/M , - - 1.8 and 
-2.7 as discussed in detail in §5.1.1 of I Vick et al.l (l2010l) . Since 
the mass loss from the surface in unseparated and the outgo- 
ing Li concentration is the Li concentration in the convection 
zone, this implies no separation in the convection zone and the 
Li mass loss flux is X(Li)pywind- It is as if the diffusion veloc- 
ity were zero there. Below the convection zone the Li flux is 
~ X(Li)p(vwind - I'd)- Element conservation at the BSCZ forces 
a strong gradient just belowrn It is seen mainly in the 0.75W1E- 
12 and 0.73W1E-12 models. Because the wind velocity hardly 
exceeds the diff'usion velocity below the bottom of the SCZ in 
the 0.68 Mq model, the gradient is less abrupt in that model. 

Finally, while VC95 obtain no effects of chemical separa- 
tion in all models with 10 '^Moyr ' (see the caption of their 
Fig. 4), lithium settling leads to ~0.15dex surface abundance 
reductions (at 12Gyr) in all our models (see Fig.|7]i. In terms 



' ' More technically, the diffusion term in vd must compensate for the 
decrease in X(Li) as one approaches the convection zone. Except where 
there is a strong X(Li) gradient, only the advective part of ^'d plays a 
significant role. Numerical accuracy required some one hundred zones 
in the small region of the strong gradient. 
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of isochrone shape, VC95 obtain a monotonous linear increase 
from the cooler stars to the hotter ones (Figs. 5 and 6 of VC95), 
while our isochrones dip in the interval 6400 K > T^ff > 6000 K 
(Fig.0. 

Some of the differences are related to the difference in mod- 
els, namely the AMbscz evolution. Other differences could in- 
clude the way the boundary condition below the SCZ was writ- 
ten and the fact that the wind term in their diffusion equation 
does not include the often overlooked flux factor (mr/M,) which 
is required in evolutionary models (see Eq. [6] of Paper I). This 
factor changes v^ind by 1 % at logAM/M, - -2 and by 10% 
at logAM/M, = -1. Given that wind and diffusion veloci- 
ties are comparable for many elements near log AM/M, = -1 
(see Figs. 3 and 4), this term modifies central concentrations. 
For instance, in the 0.779 M© model with a mass loss rate of 
5 X lO^'^Moyr"', the central hydrogen concentration decreases 
below 10~^, ~ 1.3 Gyr earlier when the correction term is in- 
cluded. The effect on surface concentrations is, however, very 
small. 



5. Comparison to models with turbulence 

Along with describing the effects of mass loss on Population II 
evolutionary models with mass loss, this paper's main objective 
is to contrast mod els with mass los s to the ones with turbulence 
presented in Rich ard et al.l (l2002bl) . In particular, one goal is to 
see if a specific mass loss rate can reproduce the surface abun- 



dances exhibited in the T6.09 mod els^ around turnoff, which are 
compatible with the Spite plateau (iRichard et al.ll2005l) . 

In Fig.|6] the surface abundances of the 0.80T6.09 and 
0.80W1E-12 models are shown to be very similar. Abundances 
of ^He and ^Li are about 0.02 dex smaller in the model with tur- 
bulence, whereas elements with Z > 12 are on average about 
0.03-0.04dex smaller (with a maximum of about 0.05 dex for 
Ar). The differences in surface abundances of elements from Be 
to Mg are smaller than the line width. 

In terms of lithium and iron isochrones, as shown in Fig.|7] 
the differences are just as small. Both types of models lead to 
nearly flat abundances (above T^s = 5800K), and their plateau 
value is nearly the same. At 12 Gyr, the maximum difference 
in Fe abundance, which occurs near Teff ^ 6100K, is about 
0.07 dex, though this could be due to the absence of turbu- 
lent models in this region. Similarly, for models hotter than 
Jeff = 5800K, the maximum difference in ^Li abundances, 
which occurs at Tf^g ^ 6200K, is 0. 1 dex, and could again 
be explained by the T6.09 isochrone's poor resolution. Below 
Tgff = 5800K, larger differences appear for lithium, likely due to 
lithium destruction in the model with turbulence. The discrep- 
ancy is greater than 0.3 dex, although again, the resolution is not 
sufficient to quantify precisely the difference. 

Whereas the surface behavior is quite similar for most 
models with T6.09 and 10 '^Moyr ', the internal solutions 
ar e a little more differe nt. By comparing Fig.|3] with Fig. 8 
of IRichard et~an (12002 b'). one notes that in the 0.80W1E-12 
model, abundances change gradually from the bottom of the 
SCZ down to about log AM/M, = -1, with most variations be- 
tween 0.1 and 0.15 dex. However, in the 0.80T6.09 model at 
turnoff, abundances are constant from the surface down to just 
below log AM/Mt - -2. Even though in the mass loss mod- 
els the abundance variations are relatively small, and the corre- 
sponding changes in opacity are likely to be small as well, this 
could potentially be tested by asteroseismic observations. 



6. Subgiant and giant branches 

The evolution of the surface abundances of Li, Fe, and Ca are 
shown in Fig.|9]for 0.779 Mq models with Zq = 0.00027B It 
is shown as a function of Tgff in order to facilitate compari- 
son with stars on the subgiant and giant branches. At turnoff, 
the T6.25 model has approximately the same Li abundance as 
the 5 X 10 '^Moyr^' model and the T6.09 model approximately 
the same as the lO^'^Moyr"' model. Past turnoff, the surface Li 
abundance first increases in both models with mass loss (with a 
maximum at T^ff ~ 5800 K) as the peak in Li abundance seen 
at log AM/M, ~ -2 in the left panel of Fig. [8] is mixed to the 
surface. In the models with turbulence, this Li peak is either re- 
duced (T6.09 m odel) or efi minated (T6.25 model) as may be 
seen in Fig. 7 of IRichard et al . (2005). In the latter model, the 

Li abundance decreases continuously past turnoff to 1.5 on 

the giant branch (Teff < 5000 K). More Li has been destroyed by 
turbulence in the T6.25 model than carried away by a mass loss 
rate of 10"'^ Moyr"' since the Li abundance on the giant branch 
is smaller in the T6.25 model than in the mass loss one. This 
implies that turbulence leads to a 0.2 dex additional destruction 
of Li in the T6.25 model when compared to models with atomic 
diffusion alone. 

These res ults may be compared to the calculations a nd ob- 
servati ons of iMucciarelli et al.l (1201 ll) and of .Nordlander et al.l 
(l2012h . The T6.25 model leads to an additional 0.2 dex destruc- 
tion of Li as appears required for the observations of Li in lower 
red giant branch stars (Mucciarelli et al., 2011) to be compati- 
ble with WMAP measurements. Neither of the mass loss models 
leads to as large a reduction of surface Li on the giant branch as 
the T6.25 model, so that the RGB observations favor the turbu- 
lent T6.25 model. However, one may also compare our results to 
the Li observations in Figs. 2 and 5 of Nordlander et al. (20121) 
where isochrones of Li abundance are compared to observations 
of NGC6397. While we do not have isochrones, a comparison 
is possible using evolutionary variations of Li abundance since 
the same turbulent models were used in our calculations and for 
their Figs. 2 and 5. If one compares the l O^'^M pyr'' results, to 
the T6.0 preferred by iNordlander et al.l (1201 2l) . one notes sim- 
ilar agreement at turnoff, but better agreement on the subgiant 
branch since maxima in our Fig. |9] appear at a slightly higher 
Teft and are more pronounced for the 10 '^Moyr ' model than 
for the two turbulent cases. The agreement is also better for the 
bRGB and RGB observations since the 10 '^Moyr ' curve is 
lower than even the T6.09 curve in our Fig.|9] It then appears 
that models with mass loss lead to at least as good agreement 
with Li observations as do models with turbulence. 

For Ca and Fe, the radiative accelerations complicate the 
evolution in the mass loss cases. For instance, the large Ca abun- 
dance close to the surface for the 1 1 .664 Gyr curve for Ca on 
Fig.|3]is caused by the grad(Ca) minimum as discussed in relation 
to that figure. In turn it appears as the abundance maximum at 
Teif ~ 6500 K in Fig.|9l The Fe local maximum at T^s ~ 6000 K 
in Fig.|9]is similarly caused by the ^iad(Fe) minimum for the 
11.664 Gyr curve at log AM/M» ~ -3 in Fig. [3] These proper- 
ties affect the comparison to observations for the mass loss mod- 
el s. The comparison is sho wn for the turbulent models in Fig. 2 
of ' Nordlander et al] (|2012|) . If one assumes the abundance shifts 
just discussed, the agreement with observations would be poorer 
with the mass loss models for both Fe and Ca at T^ff ~ 6400 K 
while at Teft ~ 6000 K, it would be poorer for Fe and unchanged 
for Ca. 



'- As is appropriate for NGC 6397. As stated in §[2]and|3] most other 
models have Zo = 0.00017. 
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Fig. 9. Evolution of the surface abundances of Li, Fe, and Ca as a function of Tetj for two models with turbulence and two with 
mass loss identified in the figure. The chosen scale emphasizes the behavior on the subgiant and bottom of the giant branch. All 
models start (pre-main-sequence) to the right of the upper horizontal line and evolve to the left until they reach the main sequence 
at Teff ~ 6100 K. Then, diffusion starts being effective and the Li abundance, for instance, is reduced by 0.2 to 0.4 dex as the star 
approaches tumoff, at Teff ~ 6600 K. The models then move to the subgiant and giant branches where dilution further reduces the 
surface Li abundance by a factor of ~ 10. For Ca and Fe the radiative acceleration complicates the surface abundance evolution in 
the mass loss cases. 



Given the difficulty justifying such a high mass loss rate we 
did not push the comparison with observations further 



7. Conclusion 

Evolutionary models with atomic diffusion and unseparated 
mass loss explain some characteristics observed in galactic halo 
dwarfs. A mass loss rate of 5 x lO"'^M0yr"' leads to surface 
lithium isochrones that exhibit plateau-like behavior with a max- 
imum scatter of about 0. 1 dex around an average value that is - 
0.3 dex below the initial lithium abundance (Fig.|7]i. For models 
with mass loss rates of lO^'^Moyr"', the lithium isochrones are 
even flatter with a scatter of about 0.05 dex around an average 
value about -0.2 dex lower than the initial abundance. 

While mass loss does not have a significant impact on a star's 
structural properties, its effects on abundances, both in the inte- 
rior and on the surface, can be important. In the interior, abun- 
dance profiles are either dominated by atomic diffusion — and re- 
flect local variations in grad — or are dominated by the wind — 
and adjust in order to conserve the flux coming from below 
(Sect. l3.11 l. A mass loss rate of 10^'^Moyr"' reduces both in- 
ternal and surface anomalies by up to a factor of 6-7, while a 
mass loss rate of 10^'^Moyr"' reduces all anomalies to below 
0.2 dex (Fig.E] and ig. While the 10 '-Moyr^' models lead to 
generalized underabundances for all elements, lower mass loss 
rates allow for overabundances to develop at the surface. 

The mass loss models prese nted in this pap er differ signif- 
icantly from those presented in IVauclair & Charbonnel (igQSh 
(Fig. IS). Although some of this difference can be attributed to 
differences in input parameters, other discrepancies are more 
difficult to explain (Sect.|4]i. In one critical aspect however, 
the results are robust: both calculations conclude that below 
lO^^'^Moyr"' mass loss has virtually no effect while above 
lO^'^Moyr"' it eliminates all effects of atomic diffusion from 
the surface. 

Wit h respect to the T6.09 models from iRichard et al.l 
( l2002bl) . the models with a mass loss rate of 10^'^Moyr"' 
are very similar above Teft = 5800K (Figs.|6] and [7]). Below 



this temperature, turbulence leads to greater lithium destruction. 
Because of differences in the internal distribution of elements 
(Sect.|5]l, asteroseismology could perhaps distinguish between 
the two scenarios. 

Unfortunately, even for mass loss rates of 10^'^Moyr"', the 
implied stellar winds are barely detectable through direct obser- 
vation. Furthermore, considering that acceleration processes for 
these types of winds are unknown, it is difficult to justify mass 
loss rates that are 50 times higher than the solar mass loss rate 
(~ 2 X 10"''*Moyr"'). It is unlikely that winds in Population II 
stars are driven by radiation; therefore, if these stars host solar- 
like winds, one could hope to observe solar-like coronae and/or 
X-ray emission. 

A signature of a chromosphere has recently been observed 
in Pop II dwarfs: the He line 10 830 A is present in almost all 
dwarfs with [Fe/H] = -0.5 to -3.7 (Takeda & Takada-Hidal 
1201 lb . Signs of chromospheres had earl ier been reported in 
these stars by Peterson & Schriiverl(l200ll) . On the other hand, 
I Cranmer & Saarl (1201 ll) have recently evaluated mass loss rates 
in various types of stars. According to their Fig. 13, mass loss 
rates of ~ 10 '^Moyr^' are to be expected in dwarf stars of 
Teff ~ 6000 K of solar metallicity that have rotation periods of 
around six days. They do not give equivalent values for lower 
metallicity stars. The most critical factor is, however, the rotation 
rate. Pop II dwarfs are believed to rotat e slowly. However, there 
are ob servations (reviewed in Fig. 6 of iTakeda & Takada-Hidail 
1201 ih suggesting that many rotate at Ve sin ; ~ 5 km/s. It would 
seem important to confirm both the rotation rate observations 
and expected mass loss rates for low metallicity stars. In the 
meantime it seems premature to disregard mass loss as a pro- 
cess competing with atomic diffusion in Pop II stars. 

As a next step, these models could be extended to the 
He flash in order to further constrain the effects of mass loss 
in Population II giants, for which mass loss has extensively 
been observed and studied (Origliaet al. 2007). It would also 
be interesting to see if such mass loss rates w ould affect, or 
even eliminate, the effects of atomic diffusion (iMichaud et al.l 
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201C ) and/or of therm ohaline mixing (ICharbonnel & ZahnI 



2007l:lTraxleretal.ll2011l) 
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Appendix A: Version of Figures|3]and|4]witli all 
atomic species included in the calculations. 
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Fig. A.l. Radiative accelerations and gravity at four different ages in a 0.8 M© model with M - 10 '^Moyr '. The vertical lines 
indicate the position of the bottom of the SCZ. The ages are identified in the figure. 
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